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Since the discovery of IGR J 1 63 1 8 — 4848 in 2003 January, INTEGRAL has accumulated more 
than 5.8 Ms in 1BIS/1SGRI. We present the first extensive analysis of the archival INTEGRAL 
data (IBIS/ISGRI, and JEM-X when available) for this source, together with the observations 
carried out by XMM-Newtoa (twice in 2003, and twice in 2004) and Suzaku (2006). 

The source is very variable in the long-tenn, with periods of low activity, where the source is al- 
most not detected, and flares with a luminosity ~ 1 0 times greater than its average value (5.4cts/s). 
IGR J 1 63 18 — 4848 is a HMXB containing a sgB[e] star and a compact object (most probably a 
neutron star) deeply embedded in the stellar wind of the mass donor. The variability of the source 
(also in the short-term) can be ascribed to the wind of the optical star being very clumpy. 

We study the variation of the spectral parameters in time scales of INTEGRAL revolutions. The 
photoelectric absorption is, with IVh around 1 0 24 cm - ’, unusually high. During brighter phases the 
strong K-alpha iron line known from XMM-Newton and Suzaku observations is also detectable 
with the JEM-X instrument. 
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Figure 1 : 29° x 9° image obtained with IBIS/ISGRI in the 40-60 keV range for the revolution 0468. 

1. IGR J16318— 4848 

IGR J 163 18 4848 was discovered on 2003 Jan 29 with the IBIS/ISGRI soft gamma-ray de- 
tector onboard INTEGRAL. It has been proposed that the source is a High Mass X-ray Binary 
(HMXB) with an sgB[e] star as the mass donor (e.g., [1, 2, 3]) surrounded by a dense and absorb- 
ing circumstellar material. The distance is estimated to be between 0.9 and 6.2 kpc. 

We analyzed the first Suzaku observation of IGR J 163 18-4848 [4], and described the average 
X-ray spectrum of the source with a strongly absorbed power law continuum with a photon index 
of 0.676(42) and an exponential cutoff at 20.5(6) keV. These parameters are typical for accreting 
neutron stars (e.g., [5, 6, 3]). We searched for pulsations in the range between 1 s and lOksec but 
could not find any coherent oscillation, so a black hole is not ruled out. Consistent with earlier 
work, strong fluorescent emission lines of Fe Ka, Fe K/3, and Ni Ka are observed, as well as a 
high column density, Mt = 1 .95(3) x 10 24 cm" 2 . 

2. Spectral Analysis 

Figure 2 shows the results of the analysis of the simultaneous data from IBIS/ISGRI and JEM- 
X, including also the four XMM-Newton observations performed during 2003 and 2004. We fitted 
the same model we used for our Suzaku observation to the INTEGRAL data. The model consists 
of an absorbed cutoff powcrlaw with three gaussians to account for the fluorescence lines, or only 
the ones for Fe Ka and K/3 in the cases where we have only JEM-X as the low-energy instrument. 
Photoabsorption was modeled with a revised version of TBabs model 1 [7], using the interstellar 
medium abundances summarized by [7]. We let the iron abundance free for the fit to Suzaku data, 
and then fixed it, as well as the energies of the lines, to fit the other observations. As the lines are 
unresolved, the width of the lines was also fixed to cr = 0.1 eV. The normalization for Fe Ka and 
K/3 was fixed to the ratio given by the fluorescence yield, Ykp/Ykcc** 0.12537 [8], 

! http://pulsar.stemwarte.uni-erlangen.de/wilms/research/tbabs/ 
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Figure 2: Spectra of all available XMM-Newton , JEM-X and Suzaku observations of IGR J 1 63 1 8—4848 
together with the simultaneous data for IBIS/ISGRI. The model fitted is an absorbed cutoff' powerlaw, with 
gaussians for the fluorescence lines (see text). 


The spectrum below 5 keV was ignored as it is most probably due to a source close to IGR 
J 1 63 18 — 4848 ([9, 10]). 

Figure 3 shows the variation of the spectral parameters. Since the parameters are better con- 
strained with XMM-Newton than with JEM-X, we will focus on the former. We can see that there 
is a change in the hydrogen column from 1.1(6) x 10 24 cm 2 to 3.4(7) x 10 24 cm 2 . We did not 
detect a Compton shoulder for the Ka line at an upper limit for its flux of 1 .8 x 1 (j 5 ph cm 2 s 1 
[4]. These two facts are in agreement with the picture of the absorbing medium being the super- 
giant’s highly inhomogeneous wind. The photon index and the cutoff energy are changing along 
the different observations as well, however, they are also consistent with a constant value within 
their uncertainties. 
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Figure 3: Parameters obtained from the fits to all available XMM-Newton, JEM-X and Suzaku observations 
together with IBIS/1SGRI (Fig. 2). The dashed lines show the (unweighted) mean of the corresponding 
values in each panel, except for the y 2 ed , where the dashed line is at 1 .0 to check for the goodness of the fit. 

3. Variability 

The light curve for the energy range 20-40 keV of IBIS-ISGR1 is shown in Fig. 4. The upper 
panel displays the total light curve with all the IB1S/ISGR1 data available to us. The lower panels 
show some close-ups on the general light curve. 

1GR J 163 18—4848 is varying significantly on short and long timescales. In the lower pan- 
els we can see how the source is constantly changing, with small short flares that are alternated 
with longer ones, like the one seen in the fourth lower panel. In these more intense flares, IGR 
J163 18—4848 reaches more than 10 times its average count rate of 5.4 cps in the 20-40 keV band 
(IB1S/1SGR1). 

This flaring behavior is typical for wind accreting systems, suggesting once more that IGR 
J 163 18-4848 is a system of a compact object imbedded in the dense wind of the sgB[e] star. 
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Figure 4: Upper panne): 1BIS/ISGRI light curve (20-40 keV) for all the archival data on IGR J 163 18-4848. 
Lower pannels: close-ups of different parts of the total light curve that show a flaring behavior typical of 
wind accreting systems (see text). 


We were aiming to detect either a pulse period or an orbital period, thus we searched for 
a period from 1 s to 900 ks. We looked for any coherent oscillation in the light curve given by 
the events detected with the XIS and PIN onboard Suzaku, and also in a 75 s-binned light curve 
generated from all archival INTEGRAL data. We used the standard epoch folding for the search, 
but the result was negative and therefore we cannot determine the nature of the compact object in 
IGR J163 18-4848. 

4. Outlook 

We will carry out a deeper analysis on the data, in particular we will: 

• look for pulsations or the orbital period in the general light curve for 1B1S/ISGRI with more 
sophisticated methods. 

• study the fluorescence lines in detail: we will look for ionization changes in the long term 
(using the 5 XMM-Newton and Suzaku observations), and in the short term, dividing these 
observations into shorter intervals. A preliminar analysis showed us that we can still con- 
strain the parameters for exposure times as short as 5 ks each. 

• analyze the IB1S/1SGR1 observations (i.e. the observations where only this instrument is 
available), defining hardness ratios to look for spectral changes and compare those to actual 
changes in the continuum parameters provided by fits. 

• compare the data from the satellites with Monte-Carlo simulations of the Compton scattering 
process. 

• extend the analysis to other IGR sources, and compare the results obtained. 
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